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Proton rapidity distributions have been measured by the NA49 collaboration in 40 and 158 
GeV/nucleon Pb+Pb collisions as function of collision centrality. We find that the shape and 
the yield per wounded nucleon in the mid-rapidity region vary little with centrality and are simi- 
lar to the distributions obtained from inelastic p+p interactions. This observation is satisfactorily 
described by the transport models HSD and UrQMD, although there are significant differences 
in the details of the spectral shape between the experimental data and the models as well as be- 
tween the models. The approximate in variance of the normalized proton spectrum in the vicinity 
of mid-rapidity suggests that multiple nucleon-nucleon interactions in nuclear collisions at SPS 
energies have little effect on the spectra of those final state protons which are slowed down the 
most. 
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1. Introduction 



It is now generally accepted that heavy ion collisions at ulti^a relativistic energies result in a 
fireball of matter with high density and temperature. Such conditions prevail, when the incoming 
nucleons deposit enough of their kinetic energy in the reaction zone. Little is known about this 
stopping process and about possible differences between the stopping of the incident nucleons in 
elementary nucleon-nucleon interactions and nucleus-nucleus collisions. On a microscopic basis 
and to first approximation each nucleon interacts only once in elementary interactions, whereas 
in central collisions between heavy nuclei it may scatter more than 4 times. The final state dis- 
tributions of nucleons will be different in the two cases. We study such distributions as function 
of centrality in Pb-i-Pb collisions in order to learn about the redistribution process of the incident 
nucleons. 

Nuclear stopping was studied in central collisions as a function of beam energy at AGS in 
Brookhaven [1], at CERN SPS in Geneva [2] and at RHIC in Brookhaven [3]. The inelasticity K 
is a convenient way to characterize the amount of stopping with a single variable. It is defined as 
K = Einei/{y/s/2 — Mp) with being equal to the average energy loss of the incident nucleons. 
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Figure 1: The inelasticity K (for the definition see 
the text) as function of ^s^n in central Pb+Pb 
(NA49), S+S (NA35), and Au+Au (BRAHMS) 
collisions [4]. Also shown is the commonly as- 
sumed value of K (=0.5) in inelastic p+p interac- 
tions (black line) and the result from a new NA49 
measurement of inclusive proton spectra in p+p 
coUisions at 158 GeV/c [5]. 
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Figure 2: Proton rapidity distributions from 
Au+Au collisions at 6, 8, and 10.8 GeV/nucleon 
for different centralities [1]. The arrows indicate 
the centers of the two Gaussians used to fit the dis- 
tributions. 



Fig. 1 shows K as function of ^/s. The inelasticity in central nuclear collisions is significantly 
larger than in p+p collisions (K=0.5) as expected [7]. In nuclear collisions it increases from 0.65 
to 0.78 when comparing AGS with RHIC data [4]. This increase goes along with a change of 
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shape of the rapidity distribution, namely a convex form at low (AGS) to a concave one at high 
(RHIC) energies. The different elasticities of small (p+p) and large (e.g. Pb+Pb) systems should 
also manifest in a different form of the final state proton rapidity distributions. Such studies have 
been already performed at AGS by the E917 collaboration at three energies [1]. Their findings are 
summarized in Fig. 2. Here the rapidity distribution evolves from a concave form in peripheral 
collisions via a double hump structure (most prominent at 8 GeV/nucleon) to a convex shape in 
central collisions. 

In this contribution we present the centrality dependence of proton rapidity distributions in 
Pb+Pb collisions at 40 and 158 GeV/nucleon which were obtained with the NA49 detector. The 
phase space coverage extends 1.6 units from mid-rapidity into the forward hemisphere and ranges 
from zero to 2.0 GeV/c in transverse momentum pj. 



2. Experimental setup and data set 

The NA49 detector is a large acceptance hadron spectrometer at the CERN SPS [9]. The main 
components are four large time projection chambers (TPCs) and two super-conducting dipole mag- 
nets with aim vertical gap, aligned in a row, and a total bending power of 9 Tm. Two 2 m long 
TPCs (VTPCs) inside the magnets allow for precise tracking, momentum determination, vertex re- 
construction and particle identification (PID) by the measurement of the energy loss {dE / dx) in the 
detector gas. The other two TPCs (MTPCs) have large dimensions (4m x 4m x 1.2m) and provide 
additional momentum resolution for high momentum particles as well as PID by dS/dx measure- 
ment with a resolution of around 4%. Two time-of-flight scintillator arrays with close to 1000 
photomultipliers each, complement particle identification in the momentum region 3-10 GeV/c 
(around mid-rapidity). A Veto Calorimeter (VCAL), which is placed further downstream along the 
beam and covers the projectile spectator phase space region, is used to select event centrality. The 
NA49-detector is described in detail in [9]. The primary Pb beam had a typical intensity of lO'* 
ions/s and impinged on a target Pb foil with a thickness of 224 mg/cm^. It passed through several 
particle detectors from which the timing signals (resolution 60 - 70 ps) and the individual beam 
particle trajectories were determined. A minimum bias trigger was derived from the signal of a gas 
Cherenkov device right behind the target. Only interactions which reduce the beam charge as seen 
by this detector by at least 10% are accepted. The interaction cross section thus defined is 5.7 barn. 
The contamination by fake peripheral events remaining after cuts on vertex position and quality 
amounts to less than 5% for the most peripheral collisions (see [10]). The resulting event ensemble 
was divided into five centrality bins CO, CI, C2, C3, and C4 (see Table 1 and [10]). The centrality 
selection is based on the forwai^d going energy of projectile spectators as measured in VCAL. The 
number of projectile spectator nucleons and the average number of interacting (wounded) nucle- 
ons (N^) were calculated from the selected cross section fractions using the Glauber approach. 
The track finding efficiency and dE/dx resolution were optimized by track quality criteria. To be 
accepted, each track must have at least 50 (out of a maximum of 90) potential points in the MT- 
PCs, and 50% of all potential points in the MTPCs must have been reconstructed. In addition, all 
accepted tracks are required to have at least 5 measured and 10 potential points in the VTPCs. 
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Table 1: Cross section fractions in %, average numbers of wounded nucleons (A^w) , numbers of ana- 
lyzed events and momentum averaged dE/dx resolutions (a) for the five centrality classes at 40 and 158 
GeV/nucleon. 



centrality bin 


centrality [ % ] 




accepted events 


(a) [ % ] 


40 GeV/nucleon 


CO 


0-5 


351 ±3 


13034 


3.6 ± 0.2 


CI 


5-12.5 


290 ±4 


22971 


3.7 ± 0.3 


C2 


12.5-23.5 


210 ±6 


34035 


3.7 ± 0.2 


C3 


23.5-33.5 


142 ± 8 


32668 


3.7 ± 0.2 


C4 


33.5-43.5 


93 ±7 


32071 


3.5 ± 0.2 


158 GeV/nucleon 


CO 


0-5 


352 ±3 


15306 


3.8 ±0.3 


CI 


5-12.5 


281 ±4 


23548 


3.9 ± 0.2 


C2 


12.5-23.5 


196 ±6 


37053 


3.8 ± 0.2 


C3 


23.5-33.5 


128 ± 8 


34554 


3.8 ± 0.4 


C4 


33.5-43.5 


85 ±7 


34583 


3.7 ± 0.3 



3. Analysis method 

The protons are identified by the measurement of their specific energy loss dS/dx in the rel- 
ativistic rise region measured in the MTPCs. Their yield was extracted by fitting the function 
F{p,pi) (see Eq. (3.1)) to the dE/dx distributions of all positively (negatively) charged particles 
in narrow bins of total momentum p and transverse momentum pj . The shape of F{p,pj) is 
assumed to be the sum of Gaussians. Their parameters depend on the particle masses and the 
measured track lengths. We modified the Gaussian functions by means of an extra asymmetry pa- 
rameter to account for tails of the Landau distributions which are still present even after truncation 
of the high dE/dx tails. 

The fit function F reads: 



i=d j:k,7i,e Lii^ii V2naij 



1 / X — X: 



2\{l±S)Oi 



ij 



(3.1) 



where 

• Ai is the raw yield of the particle / under consideration in a given phase space bin, 

• ni is the number of tracks in a given track length interval /. The second sum together with the 
normalization ^/ ni forms the weighted average of the track ensembles in each phase space 
interval, 

• f, is the expected most probable dE/dx value for particle type / and will be referred to as 
peak position. 

• X is the measured dS/dx of a single track. 
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• Oij is the width for the asymmetric Gaussian of particle type / in length interval /, and 

• 5 is the asymmetry parameter. 

The amplitude parameters Ai{p,pj) were determined by a maximum likelihood fit to the 
dE/dx distribution in each p, pj bin. The peak positions Xi{p) are considered to be pj independent 
and were determined by fits to the pj integrated distributions in each p bin and centrality class. The 
widths of the Gaussians (Oij{p)) depend on the particle type / and on the track length I according 
to Ojj = a ■ {xi/xji)^{\/ \/7). The exponent a (= 0.625) was extracted from simultaneous fits to 
distributions from TOF and to AE/Ax distributions from the TPCs [11]. The global widths o{p) 
were determined for each centrality bin by fits to the whole p range . These momentum averaged 
widths (a) turned out to be approximately 4 % for each centrality bin (see Table 1). 8 was derived 
from the observed differences between the widths of the right and left halves of the modified Gaus- 
sian functions. This parameter was fixed to be 0.071, since it did not show a significant dependence 
on centrality, total and transverse momentum. The fit parameters in each p, pj bin are the five raw 
yields A; of deuterons, protons, kaons, pions,and positrons. These yields were transformed from a 
fine grid in log{p), pj to a coarser grid in y,pT. 

3.1 Acceptance and inefficiency 

The raw particle yields have to be corrected for losses due to tracks which do not pass through 
the detectors or which do not fulfill the acceptance criteria (acceptance losses) and tracks which are 
not properly reconstructed (efficiency losses). The acceptance was calculated by generating a sam- 
ple of protons with flat distributions in transverse momentum and rapidity. The generated particles 
were propagated through the detectors (and the magnetic field) using the programs provided in the 
GEANT 3.21 [12] package. Along the resulting trajectories realistic detector signals are generated 
and processed in exactly the same way as the experimental data. The ratio of generated to accepted 
particle tracks in each y,pj bin is the acceptance correction factor. In addition to the well-defined 
losses due to the limited and constrained acceptance, the raw spectra may be subject to losses due to 
multiplicity and thus centrality dependent inefficiencies. These losses were minimized by restrict- 
ing the analysis to tracks with azimuthal emission angles within it 30 degrees with respect to the 
bending plane. The remaining losses were determined by the following procedure. Ten GEANT 
generated proton tracks and their signals were embedded into raw data of real events. Only those 
tracks ai^e embedded which pass all acceptance criteria. These modified events are reconstructed 
with the standard reconstruction chain. The ratio of all generated tracks to those reconstructed con- 
stitutes the connection factor for reconstruction inefficiencies. The resulting efficiencies vary with 
centrality by less than 5% and are over 95% in most y,pT bins. 

4. Results 



Fig. 3 presents the main experimental results of our study. Proton rapidity distributions are shown 
at 5 different centralities and two beam energies. The centrality intervals range from the - 5% 
most central events to the ensemble of events which populate the interval from 33.5 to 43.5%. We 
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Figure 3: Proton rapidity distributions from Pb+Pb collisions at 40 (left) and 158 GeV/nucleon (right) for 
different centralities.The full symbols stand for measured points.The open symbols are reflected at mid- 
rapidity. Beam rapidities are 2.2 at 40 and 2.9 at 158 GeV/nucleon. 
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Figure 4: Proton rapidity distributions from Pb+Pb collisions at 158 GeV/nucleon for different centralities. 
The full symbols stand for measured points. The open symbols are reflected at mid-rapidity. All spectra are 
normalized to the number of wounden nucleons ((A^v)). Also shown are the results of model calculations 
(HSD and UrQMD) 

note approximate shape independence at 158 GeV/nucleon (right). At 40 GeV/nucleon (left) the 
shape shows a similar evolution as at AGS (see Fig. 2) although a dip at mid-rapidity persists at all 
energies. We compare our finding at 158 GeV/nucleon with the microscopic transport models HSD 
[13] and UrQMD [14] in Fig. 4 (model calculations at 40 GeV/nucleon results were not available 
at the time of the workshop). The HSD model does remarkably well in shape and in magnitude 
at all centralities. In UrQMD the protons are slowed down too much in central collisions and too 
little in the more peripheral ones. 
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Figure 5: Proton number density as function of 
laboratory rapidity from inelastic p+p interactions 
at 158 GeV/c (red circles). The experimental data 
were measured in the forward hemisphere [5], 
but are shown here in the backward hemisphere 
(by reflection at mid-rapidity). Also shown are 
experimental results from p+p interactions at 24 
GeV/c (triangles) [6] and the average of the two 
(squares). 



Figure 6: Net neutron number density as func- 
tion of laboratory rapidity from inelastic p+p in- 
teractions at 158 GeV/c. The experimental data 
were measured in the forward hemisphere but are 
shown here in the backward hemisphere (obtained 
by reflection at mid-rapidity). 



Next we compare the proton rapidity spectra of p+p and Pb+Pb both for experimental and 
model data. We first take a look at the new NA49 results on proton distributions in inelastic p+p 
collisions at 158 GeV/c [5]. Fig. 5 shows dn/dy of net protons* obtained by NA49 as function 
of laboratory rapidity together with 24 GeV/c data [6]. The shapes of both proton distributions 
are very similar. We can thus use their average as reference for the comparison of the Pb+Pb 
data at 40 GeV/nucleon. However p+p and Pb+Pb do not have the same isospin composition. 
In fact the proper reference would be the appropriately averaged proton and neutron distributions 
from inelastic p+p, p+n, and n+n interactions. Since the latter two are not availble, we consider 
the sum of protons and neutrons from p+p as reasonable approximation. The new NA49 data on 
neutron production in p+p interactions at 158 GeV were measured as function of Xf [5]. We have 
transformed this distribution into dn/dy using the appropriate Jacobian and for the mean transverse 
masse values between 1.1 and 0.97 GeV/c^ at ycm = 0.8 and 2.8 respectively. The result is shown 
in Fig. 6. 

Before we come to the detailed comparison of net proton rapidity distributions from Pb+Pb 
and p+p colUsions we check how well the two models reproduce the new p+p data in Fig. 7. The 
overall features of the experimental data are reasonably well described by HSD and UrQMD. Sig- 

'Net protons are obtained by subtracting the anti-protons from the protons and are the appropriate observable for 
the present study. However since the the p/p ratio is small at all rapidities at 40 GeV/nucleon the subtraction has not 
been done in this dataset. We have stated explicitely "net protons" whenever subtracted spectra are presented in the 158 
GeV/nucleon dataset. 
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Figure 7: Net proton rapidity distributions from inelastic p+p interactions at 158 GeV/c [5]. The full circles 
represent the experimental data which were measured in the forward hemisphere and are reflected at mid- 
rapidity. Also shown are results of calculations with the HSD (top) [13] and UrQMD (bottom) [14] models. 



nificant differences emerge only next to the elastic peak at rapidities near to 2.5 . Here the models 
underpredict the experimental data by almost 30%. This loss is compensated by an overshoot 
around rapidities of 1.5. 

Fig. 8 summarizes the most important conclusion from the study of proton production in 
Pb+Pb and p+p collisions at SPS energies. We compare the yield of net nucleons (protons plus 
neutrons) as function of rapidity, normalised to the average number of wounded nucleons from 
p+p interactions with net protons from Pb+Pb collisions at two centralities (most central in the 
upper and most peripheral in the lower row) and two energies (40 GeV/nucleon on the left and 158 
GeV/nucleon on the right). The p+p points at the lower energy are from an interpolation of 158 
and 24 GeV/c data (see Fig. 5). We find that the Pb+Pb and p+p spectra are remarkably similar in 
shape and in yield. In particular- there is almost perfect agreement between Pb+Pb and p+p data 
in peripheral collisions at the higher energy (lower right). In central collisions the parabolic shape 
gets more shallow but the integral of the normalized measured spectrum remains almost unchanged 
(0.39 in central and 0.41 in peripheral collisions). We note here that this integral over the complete 
rapidity interval would be two, if the net baryon content in the hyperons and in bound nucleons 
is neglected. It seems that the effects due to secondary interactions of the participant nucleons 
happen mostly outside of the experimental acceptance at the higher energy. The multiple collision 
effects are coming out clearer at the lower energy (left column). At 40 GeV/nucleon a rapid, sig- 
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Figure 8: Normalized net proton rapidity distributions from Pb+Pb collisions at 40 (left) and 158 
GeV/nucleon (right) for the most central (upper row) and the most peripheral centrality selection (lower 
row). The normalization constant is (A^w)/2 to account for the unseen neutrons in the Pb+Pb data. Also 
shown as reference are the rapidity distributions of net protons plus neutrons from inelastic p+p interactions 
(open symbols) . 



nificant change in shape occurs roughly one unit away from mid-rapidity or 1.2 units away from 
target-rapidity. We note that this region (1.2 units away from target-rapidity) is not covered by our 
acceptance at the higher energy (158 GeV/nucIeon). The effect observed at 40 GeV/nucIeon could 
thus also be present at 158 GeV/nucleon. The integrated yields of the normalized 40 GeV/nucleon 
spectra are again remarkably constant (0.36 in central and 0.38 in peripheral collisions). The com- 
parison to our interpolated p+p reference is less conclusive than at 158 GeV/nucleon. It is certainly 
desirable to get the corresponding experimental data from the forthcoming NA61 experiment (see 
contribution by G. Stefanek to these proceedings). 
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In view of the significant differences in the number of colUsions per participating nucleon 
between p+p and Pb+Pb as well as in peripheral and central (Pb+Pb) collisions our findings suggest 
that multiple collisions have little influence on the protons close to mid-rapidity, i.e. those which 
experience strong stopping. This finding is more pronounced at 158 than at 40 GeV/nucleon. It 
could be verified by a detailed study of transport model caculations which seem to reproduce these 
phenomena. 
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